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Abstract 

We consider a wideband spectrum sharing system where a secondary user can share a number of orthogonal 
frequency bands where each band is licensed to an individual primary user. We address the problem of optimum 
secondary transmit power allocation for its ergodic capacity maximization subject to an average sum (across the 
bands) transmit power constraint and individual average interference constraints on the primary users. The major 
contribution of our work lies in considering quantized channel state information (CSl)(for the vector channel 
space consisting of all secondary-to-secondary and secondary-to-primary channels) at the secondary transmitter as 
opposed to the prevalent assumption of full CSI in most existing work. It is assumed that a band manager or a 
cognitive radio service provider has access to the full CSI information from the secondary and primary receivers 
and designs (offline) an optimal power codebook based on the statistical information (channel distributions) of the 
channels and feeds back the index of the codebook to the secondary transmitter for every channel realization in real- 
time, via a delay-free noiseless limited feedback channel. A modified Generalized Lloyds-type algorithm (GLA) 
is designed for deriving the optimal power codebook, which is proved to be globally convergent and empirically 
consistent. An approximate quantized power allocation (AQPA) algorithm is also presented, that performs very 
close to its GLA based counterpart for large number of feedback bits and is significantly faster. We also present 
an extension of the modified GLA based quantized power codebook design algorithm for the case when the 
feedback channel is noisy. Numerical studies illustrate that with only 3-4 bits of feedback, the modified GLA 
based algorithms provide secondary ergodic capacity very close to that achieved by full CSI and with only as little 
as 4 bits of feedback, AQPA provides a comparable performance, thus making it an attractive choice for practical 
implementation. 



I. Introduction 

Radio spectrum is a limited and precious natural resource, which, traditionally, is licensed to users by 
regulatory authorities in a very rigid manner where in order to avoid interference, the licensed owner has 
an exclusive right to access the allocated frequency band [Q]|. Consequently, as the number of wireless 
communication systems and services grows, the availability of vacant spectrum becomes severely scarce. 
However, recent measurements by the Federal Communications Commission reveal that many portions 
of spectrum are mostly under utilized or even unoccupied. This led to the idea of cognitive radio (CR) 
technology, originally introduced by J. Mitola [2], which holds tremendous promise to dramatically 
improve the efficiency of spectral utilization. The key idea behind CR is that an unlicensed/secondary 
user (SU) is allowed to communicate over the frequency band originally licensed to a primary user 
(PU), as long as the transmission of SU does not generate unfavorable impact on the operation of PU. 
Effectively, three categories of CR network paradigms have been proposed: interweave, overlay, and 
underlay [3J. In the underlay systems, which is the focus of this paper, the SU can transmit even when 
the PU is present, but the transmitted power of SU should be controlled properly so as to ensure that 
the resulting interference does not degrade the received signal quality of PU to an undesirable level (6) 
by imposing the so called interference temperature [1] constraints at PU (average or peak interference 
power (AIP/PIP) constraint). This type of CR is also known as the 'spectrum sharing' AH model. 

[|7l first studied the behavior of capacities of different AWGN channels under received-power constraints 
(AIP) at the PU receiver (PU-RX), which showed that for point-to-point non-fading AWGN channels, 
the capacity performance with transmit and received power constraints are very similar. The ergodic 
capacity of narrow band spectrum sharing model with one SU and one or multiple PU under either 
AIP or PIP constraint at PU-RX in various fading environments was studied by d), illustrating that in 
a fading environment, spectrum access opportunity for the SU significantly increases compared to the 
AWGN case. In [9], the authors studied optimum power allocation for three different capacity notions 
under both AIP and PIP constraints. [6] designed optimal power transmission strategies for maximizing 
ergodic capacity and outage capacity under various combinations of secondary transmit power constraints 
and interference constraints. 

Most of the above results assume perfect knowledge of full channel state information (CSI) including 
the SU-TX to PU-RX channels, which is hard to realize in practice. A few recent papers have emerged 
that address this concern by investigating capacity analysis with imperfect CSI. The effect of imperfect 



2 



channel estimation in the secondary to primary channels has been investigated in IfTOl by considering the 
channel estimate as a noisy version of the true CSI, and ll20l proposed a practical design paradigm for 
cognitive beamforming based on finite-rate cooperative feedback from the PU-RX to the SU-TX. Another 
recent work lfT2l also considers imperfect CSI for the SU-TX to PU-RX channel in the form of noisy 
channel estimate and quantized channel information and investigates the effect of such imperfect CSI on 
the capacity performance of the secondary user, while assuming that the SU-Tx has full knowledge of the 
SU-Tx to SU-Rx channel. Finally, ifTTll studies the issue of channel quantization for resource allocation via 
the framework of utility maximization in OFDMA based cognitive radio networks, but does not investigate 
the joint channel partitioning and rate/power codebook design problem. Indeed, the lack of a rigorous and 
systematic design methodology for quantized resource allocation algorithms in the context of cognitive 
radio networks forms the key motivation for our work. In this paper, we investigate an ergodic capacity 
optimization problem for the secondary user where quantized information about the vector channel space 
consisting of SU-Tx to SU-Rx channels and SU-Tx to PU-Rx channels is available to the SU-Tx via a 
limited feedback channel without delay. We consider a wideband spectrum sharing system where one SU 
shares M different frequency bands with M PU's, each PU using a separate band. We address the problem 
of ergodic capacity maximization of the secondary user subject to an average sum (across the bands) 
transmit power constraint on the secondary user and individual average interference constraints on the 
primary users, using quantized channel information. To this end, we assume the availability of an entity 
called a band manager (or a CR service provider) who has access to the full CSI including all secondary- 
to- secondary and secondary-to-primary channels. It designs (offline) an optimal power codebook based 
on the statistical information (channel distributions) of the channels and in real-time, feeds back the index 
of the codebook to the secondary transmitter for every channel realization, via a limited feedback link. 
The secondary transmitter then uses the corresponding power code vector for its transmission. 

We make the following key contributions: (1) We first present, very briefly, a systematic algorithm 
for optimal power allocation with full channel side information (CSI) at the secondary transmitter. This 
is a minor extension of the results in [[6l to the multiple PU case. However, the novelty lies in exactly 
characterizing the optimal power allocation policy based on the relationship between the available total 
average SU transmit power and the individual average interference levels at the PU receivers. (2) Next, we 
present a modified Generalized Lloyd's type algorithm (GLA) for designing the optimal power codebook 
using quantized channel information. For easier exposition, we focus on the narrowband case first and 
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present the quantized power allocation algorithm, where we prove that the modified GLA based power 
codebook design algorithm is globally convergent and empirically consistent. We provide a number 
of useful and interesting properties of the quantized powers. Then we present a complete description 
of the optimal power codebook design algorithm for the wideband spectrum sharing case under the 
average transmit power and average interference constraints. We believe this paper is the first to provide 
a systematic quantized power allocation algorithm with limited feedback for the spectrum sharing scenario 
in cognitive radio. (3) Although an offline algorithm, GLA based quantizer designs usually require a large 
number of training samples and can be computationally expensive. We therefore design an approximate 
quantized power allocation algorithm based on the derived properties of the power codebook, which is 
computationally much faster. (4) We then generalize the modified GLA based algorithm for quantized 
power allocation algorithm to the case where the limited feedback channel is noisy but memory less. (5) 
We present a comprehensive set of numerical results that illustrate (i) how the modified GLA-based power 
codebook can achieve a secondary ergodic capacity with only 3-4 bits of feedback, that is very close 
to the capacity with full CSI, (ii) how the performance of the approximate quantized power allocation 
algorithm is almost indistinguishable from that of the GLA-based algorithm with B > 4 bits of feedback 
and (ii) how the performance of the quantized power allocation degrades when the noisy feedback channel 
error probability increases. 

The rest of the paper is organized as follows. Section III] presents the system model and assumptions 



about the spectrum sharing problem with limited feedback. In Section III we present the optimal power 
allocation policy when the secondary transmitter has full CSI and discuss various special cases. In Section 



IV we present the modified GLA based quantized power codebook design algorithms for the narrowband 



case followed by the wideband case. We present results on global convergence and empirical consistency 
of the GLA based algorithms and some prove some useful properties of the quantized power code vectors. 
These properties are then used to design an approximate quantized power allocation algorithm suitable 
for moderate to large number of feedback bits that has a much faster convergence time compared to its 
GLA counterpart. In Section |V} we provide a modified GLA based power codebook design algorithm 



for a noisy limited feedback channel model. Numerical results are presented in Section VI and finally, 



concluding remarks and possible extensions are presented in Section |VII[ All proofs are relegated to the 
Appendix unless otherwise mentioned. 
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II. System Model and Problem Formulation 

We consider a wideband spectrum sharing scenario with one SU and Multiple PUs, as shown in Fig. 
[TJ where a SU is allowed to use M parallel orthogonal frequency bands (Bandi to Bandu) which are 
individually licensed to PU%, . . . , PUm respectively. Regardless of the on/off status of PUi, SU uses 
the 2-th channel as long as the impact of the secondary transmission does not substantially degrade the 
received signal quality PUi. It is assumed that the the channels between the secondary transmitter (SU- 
TX) and secondary (SU-RX) receiver and those between the secondary transmitter and the each primary 
receiver are all block fading additive white Gaussian noise (BF-AWGN) channels. Let g l E H + and 
g\ E IR + denote the real-valued instantaneous channel power gains for the link between the SU-TX and 
the receiver of PUi and i-th channel between the SU-TX and SU-RX, respectively, where IR + denotes 
the set of nonnegative real numbers. These channels are assumed to be stationary ergodic with absolutely 
continuous probability density functions (pdf) /o(fl'o) an d h{9x)- F° r analytical simplicity, the interference 
from PUi-TX to SU-RX is neglected (similarly as in [[TJ, [0). In the case where the interference caused 
by the primary transmitter at the secondary receiver is significant, the SU ergodic capacity results derived 
in this paper can be taken as upper bounds on the actual capacity under primary-induced interference. All 
g l and g\ (i = 1, . . . , M) are statistically mutually independent and, without loss of generality {w.l.o.g), 
are assumed to have unity mean. Similarly, additive noises for each channel are independent Gaussian 
random variables with zero mean and unit variance w.l.o.g. When M = 1, this system becomes a typical 
narrowband spectrum sharing model considered in Q]|[[5]|[[6]|. 

Given a channel realization g = {g^, . . . ,g^} and g 1 — {g\, . . . ,gf f }, we assume that a channel 
side information (CSI) r/(go, gi) is available at the SU-TX. The power allocated at the SU-TX on the 
M parallel SU links is represented by the vector p(ry(g ,g 1 )) = {pi(^(g ,gi)), • • • ,Pm(v(So^i))}, the 
ergodic capacity of the SU for this wideband spectrum sharing system can be expressed as 

1 - 

C =mY, E ^+9\pM^Zx)))\ (1) 

i=l 

where, for simplicity, we have ignored the factor | at the front of the capacity expression and log 
represents the natural logarithm. A common way to protect PU's received signal quality is by imposing 
either an average or a peak interference power (AIP/PIP) constraint at PU-RX flTJ [[5ll [J6J , although other 
forms of PU quality of service constraint such as PU's capacity loss and PU's outage probability ll2Tfl . 
It was shown in [5] that an AIP constraint is more favorable than a peak constraint especially in the 
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context of transmission over fading channels, since the AIP constraint is more flexible and can achieve 
larger SU capacity results with less PU capacity loss than those achieved by PIP. 

Motivated by this observation, we consider the following optimal power allocation scheme that maxi- 
mizes the ergodic capacity of SU in a wideband spectrum sharing scenario, under an AIP constraint at 
each Pf/j-RX and an average sum transmit power constraint (ATP) for the SU, given by, 

M 



lliax Tfj^Mi + tfiftMfo.gi)))] 



i=l 

M 



s.t. E[g l oPi ( V (g , gl ))] < Ql vg , Vz, ^ E \pM8o, Si))] < P«,g (2) 

i=i 

In the next section, we present the optimal power allocation results assuming that full channel state 
information (CSI) is available at the SU-Tx (i.e, ^(g ,gi) = (g ?gi)X followed by the case of quantized 
channel information (or limited feedback) in Section IV, where ^(g ,gi) represents a deterministic index 
mapping scheme, such that r/(g , g-J = j, j e {1,2,..., L} when the instantaneous channel gains (g , g-J 
belong to a carefully constructed partition TZj of the channel space IR^ x IR+ f . 

III. Optimal Power Allocation with Perfect Channel State Information 
In this section, we assume that SU-TX has perfect knowledge of g and g 1 (full CSI at the transmitter), 
that is, ?7(go>gi) = (go>8i)- ft i s easv t0 verify that the problem given in ^ is a convex optimization 
problem. By applying the necessary and sufficient Karush-Kuhn-Tucker (KKT) conditions for optimality, 
the optimal power allocation can be easily shown to be 



P;(g0,gl) = ( TT- J i (3) 

A/ + v\gl 9i 

where and fi{ are the nonnegative Lagrange multipliers associated with the ATP constraint and the AIP 
constraint of PUi respectively, and (x) + = max(x,0). This solution is clearly a minor extension of the 
narrowband result in [6J. However, in the wideband case (M > 1), it should be noted that determining 
the optimal power allocation scheme involves obtaining the optimal values of the (M + 1) Lagrange 
multipliers. Since all the constraints in Problem ([2]) may not hold with equality simultaneously, it is 
difficult to determine A^" and /if, Vi. Although they can be obtained by, e.g., the ellipsoid method IfTTl . 
this procedure can be time consuming. Thus motivated, we present a complete solution to Problem ([2]), 
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summarized in the following theorem (here the term "iff refers to "if and only if) 

Theorem 1: With perfect channel information ^(gf^g-J = (g^gx) at the SU-TX, the optimal power 
allocation for problem ([2]) is given by 



Pi (80.81, 



1 



( 

(± - 1 



1 \ + 



9\ 



iff Eft 



j_ 1 



Si' 



otherwise 



iff Pavg > m Y^i=l -^[(^j g\. 



<Q 



avg 



otherwise 



Proof: See Appendix A for a proof. 
One can also easily obtain the following special cases (we do not provide the proofs due to space 
constraints as they are straightforward): 



1) When M = l(narrowband spectrum sharing case), from theorem 1 the condition E[(jj 



< 



Qavg) becomes P avg < Q avg (note that we have removed the superscript from Q avg as there is only 
one primary user), and the optimal power allocation solution specialises to the one presented in 

0. 



2) When ^ = Vi, 

P*(go,8i) = (t*7 - \) + , - \) + ] < Q\ vg V. (4) 

where is given by ~ ^~) + ] = Pavg- For this case, if additionally g x are inde- 

pendent and identically distributed, we can simplify the condition E[(-k- — \) + ] < Q l avg Vi as 

Pavg < m i n {Qavgi ' ' • ) Qavg)- 

3) If Q l avg = ■ ■ ■ = Qal g = Qavg and both g and g x are independent and identically distributed, the 
optimal power allocation policy is to assign equal power to each SU link, which is identical to the 
power allocation policy for the M — 1 case. 

Appealing to the convexity of Problem ([2]), one can show that in Theorem [TJ one of the cases must 
hold, and the corresponding power allocation scheme must be the global optimal solution for the original 
problem (S. An algorithm can be then easily designed to obtain p*(g , gx), and the associated non-zero 



Lagrange multipliers can be obtained by solving the KKT optimality conditions numerically (e.g, via a 
bisection search). 

IV. Optimum Quantized Power Control with Finite-Rate Feedback 

The assumption of full CSI at the SU-TX (especially that of g ) is usually unrealistic in practical 
systems. In this section, we are therefore interested in designing power allocation schemes based on 
quantized (g , g x ) information acquired via a no-delay and error-free feedback link with limited rate. 
Here we assume that there is an entity (such as CR service provider or a band manager [4J) who can 
obtain perfect information on g 1 from SU-RX or SU base stations and perfect information on g from 
PU base stations, presumably over a wired link, and then forward some appropriately quantized CSI 
to SU-TX (and SU-RX for decoding purposes) through the feedback link. More specifically, given B 
bits of feedback, a power codebook V= {Pi, . . . , Pl} (where Pj = {pij, . . . ,pMj}, j = 1, • • • , L) of 
cardinality L = 2 B , is designed off line purely on the basis of the statistics of g , g v This codebook 
is known a priori by both SU-TX and SU-RX. The vector space of (g , g x ), is thus partitioned into L 
regions IZi, . . . ,1Zl using a quantizer Q (codebook element P^ represents the power level used in IZj 
). The CR service provider/band manager maps the current instantaneous (g , g : ) information into one 
of L integer indices and sends the corresponding index to the SU-TX via the feedback link (e.g., if the 
current (g , g x ) falls in IZj, then r/(g , g x ) = j will be conveyed back to SU-TX). The SU-TX will use 
the associated power codebook element (e.g., if the feedback signal is j, then P., will be used as the 
transmission power) to adapt its transmission strategy. 

Let Pr(TZj), E[»\lZj] denote Pr((g ,g 1 ) e TZj) (the probability that (g ,gi) falls in the region TZj) 
and -£?[•! (go, gi) £ T^j], respectively. Then the secondary ergodic capacity maximization problem ^ with 
limited feedback can be formulated as 

L M 

J - ' 3 j=l i=l 

L L M 

s.t. ^^PMPr^) < Ql g , Vz, J>,)Pr(7^) < P avg (5) 

j=\ j=i i=i 

Our objective is thus the joint optimization of the channel partition regions and the power codebook 
such that the ergodic capacity of SU is maximized under the above average transmit power and average 
interference constrains. 
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A. Narrowband spectrum-sharing case 

For ease of exposition, we first look at the relatively simpler case of M = 1 (where SU shares a 
narrowband spectrum with only one PU). For simplicity (with some abuse of notation), let pj, gi,go, Qavg 
represent p\j, g{, g\, Q\ vg respectively. Thus problem (|5j) with M — 1 becomes, 

L 

max. ^ J E[log(l + W )|^]Pr(^) 

L L 

s.t. ^ElgoPjlKjlPrCRj) < Q avg , ^ ft Pr(7^) < P avg (6) 

3=1 j=l 

We solve the problem ([6]) based on the Lagrange duality method. First we write the Lagrangian of above 
problem as 

L 

L(P,X,p) =J2 E i lo &( 1 + 9iP j ) - Xpj 
i=i 

where A and p are the nonnegative Lagrange multipliers associated with the ATP constraint and AIP 
constraint respectively. The Lagrange dual function g(X, p) is defined as 

L 

max ^[log(l + gxpj) - Xpj - pg pj |^]Pr(^) (8) 

and the corresponding dual problem is minA>o, ^>og(X,p) + XP avg + pQavg- 

We first consider solving the optimization problem ([8]) with fixed A and fi. To this end, we employ 
an algorithm similar to a Generalized Lloyd Algorithm (GLA) |[T3l . Ifl4ll to design an optimal codebook 
for problem ([8]), which is based on two optimality conditions : 1) optimum channel partitioning for a 
given codebook, also called the nearest neighbor condition (NNC) in the context of traditional vector 
quantization (VQ), and 2) optimum codebook design for a given partition, also known as the centroid 
condition (CC) (in the context of VQ) [TT4ll . GLA is usually initialized with a random choice of codebook, 
and then the above two conditions are iterated until some pre-specified convergence criterion is met. The 
same procedure is used here for designing an optimal quantizer Q, but the design criterion for our case is 
minimizing the difference between the capacity with perfect CSI and the capacity with quantized power 
allocation under the given constraints. This amounts to designing an optimal power codebook Q that 
maximizes the Lagrangian function for quantized CSI, Yl!j=i ^P-°g(l + 9iPj) ~ ^Pj ~ ^9oPj\P-j\P r {T^'j)- 
We call the corresponding quantized power allocation algorithm for a given X, p as a modified GLA. 
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In practice, this modified GLA is implemented using a sufficiently large number of training samples 
(channel realizations for go,9i)- Beginning with a random initial codebook, one can design the optimal 
partitions using the fact that the optimal partitions satisfy TZj = {(go, gi) ■ (log(l + giPj) — Xpj — l^goPj) > 
(log(l + giPn) — Xp n — /i(?oPn), ^ j} where TZj is the corresponding partition region for power level 
Pj in the codebook, and ties are broken arbitrarily. Once the optimal partitions are designed, the new 
optimal power codebook is found by solving for argmax pj . >0 E[log(l + giPj) — Xpj — fig pj\TZj]Pr(TZj), 
Vj = 1, 2, . . . , L. Given a partition, this optimization problem is convex and by using the KKT conditions, 
one can obtain the optimal power as max(j»*,0), where p* is the solution to the equation -5'[ 1+ g 1 ip . — 
(A + fig )\7Zj] = 0. These two steps are repeated until the resulting ergodic capacity converges within a 
prespecified accuracy. One needs to note that GLA cannot in general guarantee global optimality, since the 
two optimality conditions (NNC and CC) mentioned above are just necessary conditions 031. Thus it is 
possible that the our resulting quantizer is only locally optimal. While convergence of our modified GLA 
follows immediately by noting that the Lagrangian Y^=i E[log(l+giPj) — Xpj— figoPj\TZj]Pr(TZj) is non- 
decreasing at each iteration and is upper bounded (due to the finite average transmit power and average 
interference constraints), it is important and instructive to state a more formal result along the lines of lPT5l . 
Since GLA is initialized with a random codebook and the optimal partitions and codevectors are found 
using training samples drawn from empirical distributions, it is crucial that GLA is globally convergent 
with respect to the choice of initial codebooks and empirically consistent. For more formal definitions 
of these two properties, see HT31I . Under the assumption of absolutely continuous fading distributions for 
go,gi and mild regularity assumptions satisfied by these distributions, one can show that the modified 
GLA satisfies the conditions for global convergence and empirical consistency stated in 0"5l and thus we 
have the following result: 

Theorem 2: The modified GLA that solves the optimization problem ([8]) satisfies the global convergence 
and empirical consistency properties of [fT5l . 

Proof: See Appendix A for a proof of this result. 

Next, we present some useful properties of the optimal power solutions obtained via the modified GLA. 
We use the partitions TZ X , . . . , TZ L and the corresponding power levels p x , . . . ,pi to denote the convergent 
optimal solutions. 

Lemma 1: Given partitions 7£i, . . . , TZl and the corresponding power level pi, . . . ,Pl, (where TZj and 



10 



lZj + i,Vj E {1, . . . , L — 1} are adjacent regions and Pj ^ Pj+i), the boundary between any two adjacent 
regions TZj and is given by, g\ = — 6 + ^t>-+lao)(Tj-'p^ + i) which, when fi 7^ 0, is a monotonically 

log(^L) 

increasing convex function of go and as g\ — V 00, go — > ~( p -p^ — 

Proof: From the NNC condition of the modified GLA, the boundary between two adjacent regions 
TZj and TZ j+x satisfies log(l + g x pj) - Xpj - fig pj = log(l + giPj+i) - \p j+ i - figoPj+i- Solving the 
above equation for gi, the result in the above Lemma follows. It is straightforward to show that it is an 
increasing convex function of g by investigating the first and second derivatives. 



Remark 1: In case A > 0,fi = 0, the AIP constraint is inactive and the ATP constraint is satisfied 
with equality. In this case, the boundary between any two adjacent regions TZj and Hj+i becomes 

gi = — x(p - P ) • Clearly, Problem (5) reduces to an ergodic capacity maximization problem 

with quantized channel information. For the narrowband case, it becomes a scalar quantization problem 
involving quantizing g\ only. Note that while for the narrowband case, this no longer pertains to a cognitive 
radio problem, the properties of the optimal quantized power allocation scheme are still important for 
the wideband case (M > 1). This is due to the fact that in the wideband case, it is possible that for a 



specific (say the i-th) channel, the AIP constraint is inactive (p,i > 0) while A > 0. See Section |TV-B| for 
further details. 

We now give an example to illustrate what the optimum partition regions actually look like. For this 
example, g and g\ are both exponentially distributed (Rayleigh fading) with unit mean and L = 4 (2 
bits of feedback). The optimum partition regions are as shown in Fig. [2] for A > 0, /x > 0, and Fig. [3] for 

X> 0,fi = 0. 

We obtain the following properties for the optimal quantized power levels where (as illustrated in 
Figure [2]) the regions 7Zi,7Z 2 , ■ ■ ■ etc. are sequentially numbered, with 1Z\ being the region closest to the 
<7i axis and TZl being the region closest to the go axis. Note that these properties apply regardless of 
whether [i > or fi = 0. 

Theorem 3: i). P\> • • • > Pl 

ii) . All boundaries between any two adjacent partitions satisfy g\ > A + figo- 

iii) . Given B bits of feedback (or L = 2 B regions), for the first L-l regions, we always have strictly 

positive power, i.e. pi > ■ ■ ■ > Pl-x > 0, whereas pi is simply nonnegative, i.e. p L > 0. 

iv) . When A + fi > 1 (note that if A = 0, fi > 1 implies Q avg < 1, and if /i = 0, A > 1 corresponds 
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to P avg < 1), we always have pi = 0. In addition, when L (the number of quantized regions) is 
sufficiently large, no matter what A, /i is, px must be 0. Additionally, as L — > oo the boundary 
between 1Zl-i and 1Z L approaches g x = \ + /j,g and lim^oo p L -i = 0. 
Proof: See Proof in Appendices B-E. 



Remark 2: The above properties of optimal quantized power values are interesting for two reasons. 
From property ii), it is clear that (g , g t ) ElZj for j = 1, 2, . . . , L — 1 satisfy the property g x > A + fj,g 
whereas for the region IZl, this property may or may not be satisfied. Since the quantized power values in 
the first L — 1 regions are strictly positive, it is easy to relate this property to the corresponding property 
of the full CSI based optimal power value which is strictly positive if and only if when g\ > \f + fi^go. 
Also, as L — > oo, the boundary between 1Zl-i and IZl approaches g\ = A + fxg , thus making this 
relationship between the quantized power allocation scheme and the full CSI power allocation scheme 
stronger. 

Finally, property iv) allows one to obtain an approximate quantized power allocation scheme (AQPA) 
for large L by setting p L = and taking the limit as Pl~i — > 0. This is particularly useful as the modified 
GLA becomes computationally intensive for large L, whereas AQPA provides a performance that is 
extremely close to that of the modified GLA, while requiring very little computation time. A detailed 



description of the AQPA is provided in Section |IV-C| followed by illustrative numerical simulations in 
Section I 



Based on the above Lemmas, one can solve for the optimal quantized power values given a partition 
IZi, IZ2, ■ ■ ■ ,7Zl is equivalent to solving the following set of nonlinear equations for pi,P2,P3, ■ ■ ■ ,Pl- 

£[— ^ (A + =0, j = 1, ■ ■ ■ , L, p L = max(0,p i ) (9) 

-L 1 9iPj 

where if /* ^ 0, E[^ - (A + m)Wi = fr-Ji^ 7 ( A + ^o))f(9o)f(9i)dg d gi , with c 3 = 

Pj—pj + ie ^3 ^j+l> J M Pj Pj+l 

When n = 0, E[j^- - (A + = J^ij^p- ~ *)f(9i)dgi, with c = 00. g can be solved 

efficiently by any suitable nonlinear equation solver. 

Now that we have an algorithm based on the modified GLA for solving for the (possibly locally 
optimal) quantized power values for fixed A,/i, we can go back to solving the dual problem for finding 
the optimal values A and /i. To this end, we solve the associated KKT conditions (involving the average 
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power and the average interference constraints) numerically (e.g, via a bisection method). One can thus 
repeat the above two steps by solving ([8]) and the dual problem iteratively until a satisfactory convergence 
criterion is met. An algorithmic format for this procedure is provided for the more general wideband 
(M > 1) case in the next subsection. 



B. Wideband spectrum-sharing case 

The above algorithm for the narrowband case can be easily extended to the wideband case corresponding 
to the original problem ([2]). For this scenario, the Lagrangian function becomes, 

L AI 

L(P, A, u) = £( — £ £[log(l + gipij)\Kj])Pr(K3) 

3=1 »=1 

L M ML 

3=1 i—1 i—1 j — 1 

where A and [L{ are the nonnegative Lagrange multipliers associated with the ATP constraint and ith AIP 
constraint respectively. The Lagrange dual function g(\, {/^}) is defined as 

M L 



max hll E t l0 S( 1 + - - faoPiilKAPriKj) ( n ) 



i=i j=i 



where /4 = Mfa, Vi, and the dual problem is min A > 0j ,4>o,vi0(A, {//■}) + AP ai , 3 + £\ =1 ffQa 



aug- 



Similar to the narrowband case, we first consider the problem (11) to obtain g(X, with given 

A and Denote by 7^*- the j-th quantization region for the i-th band where |jf=i = ^i- Then 



problem ( 11 ) can be decomposed into M parallel subproblems, where for each band i, % = 1, . . . , M 

L 

max V £[log(l + g\ Pij ) - X Pij - ^gl Plj \n)]Pr{n)) (12) 

3=1 

is defined as the sub-dual function g^X,^) and g(X,{bt' i }) = -^Y^!=i9i{\ fA)- This kind of duality 
method is also known as the 'dual decomposition algorithm' [fT6ll . Since each subproblem ( [12] ) is similar 
to the problem ([8]) for the narrowband case and can be similarly solved by using a modified GLA. A 
and {/jfj} can be also obtained in a manner similar to the narrowband case. These two steps are then 
repeated until a satisfactory convergence criterion is met. Due to the increased complexity resulting from 
the presence of multiple bands, we provide below a description of the overall optimization algorithm 
(Algorithm 1) for solving Q. 
Algorithm 1: 
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1) Let A = 0, then all p[,i = 1,...,M must satisfy ^ > 0. Starting with some random initial 
power codebook, for each i, find /i- by solving J2j=i ^[9oPij\^j]P r C^j) = Qavg an ^ then obtain 
the corresponding (locally) optimal power codebook {pa, . . . ,pn,} using a modified GLA. Repeat 
these two steps until convergence resulting in a power codebook {P 1; . . . , P^}. With this codebook, 
if J2j=i(~M ^t=i ^[Piji^jD-^M^j) — Pavg, it is an optimal power codebook and stop; otherwise 
go to step 2). 

2) If 1) is not satisfied, we must have A > 0. For a given A, for each i, use the modified GLA to find 
an optimal power codebook first with //• = 0. If J2j=i ^[9oPij\^j}P r (^-j) < Qavg> then the cor- 
responding optimal codebook {p a , . . . ,pu} (obtained via the modified GLA) is an optimal solution 
for this z-th subproblem, otherwise, fi[ > 0, and can be found by solving J2j=i E[g} > Pij\'R'j\Pr(llj) = 
Q % avg . Find the corresponding optimal codebook entry {pn, . . . ,Pil} for the i-th subband , and then 
use this codebook to find an updated value of A by solving X^=i(if ^\Pij\^j})P r {^-j) = 
P avg . Repeat these steps until convergence and the final codebook will be an optimal codebook for 
the wideband spectrum sharing problem (|5]). 

Remark 3: Note that it is straightforward to extend the global convergence and empirical consistency 
results of Theorem [2] to the wideband case. Similarly, Lemma [3] also holds for the wideband case in 
the sense that the properties i)-iv) hold for each {pu,Pi2, ■ ■ -Pil}, Vz = 1, 2, . . . , M with fj, replaced by 
= 1,2, ... ,M and A representing the Lagrange multiplier associated with the average sum power 



constraint in (10). 



C. Approximate Quantized Power Allocation Algorithm (AQPA) 

Although an offline algorithm, the complexity of modified GLA for determining the optimal quantized 
power is very high for even a moderately large value of L. This is due to the fact that the optimal channel 
partitions and the corresponding optimal power codebook are obtained via empirically generating a large 
number of channel realizations as training samples. As L increases, the number of training samples 
required will also increase. Thus motivated, we use part iv) of Lemma [3] to derive a low-complexity 
suboptimal scheme for implementing the modified GLA for large L values. Below we describe this 
scheme for the narrowband case. A similar scheme for the wideband case can be designed accordingly. 

Note that part iv) of Lemma [3] states that as L — > oo, pl = and pl-i — > 0. Applying these 
approximations to Q allows us to obtain an approximate but computationally efficient algorithm (called 
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approximate quantized power allocation algorithm (AQPA)) for large L. AQPA first solves E[ j^-^ - — 
(A + fj,g )\TlL-i] = for pi-2 by substituting p L = and taking the limit Pl-i — > 0, which, if \i > 0, is 

91 -A 

equivalent to solving J A °° / i og(1+glPL _ 2 ) (ffi - (A + ngo))f(go)f(gi)dg dgi = for p L _ 2 - When // = 0, 

e ^L-2_! 

it is equivalent to solving for Pl-2 from J A PL ~ 2 (grj — \)f(gi)dgi = 0. Note that the above equations (for 
both n > and fi = 0) involve only one variable: £>l_2 and are thus straightforward to solve. One can then 
recursively compute £>l-3,Pl-4, • • • , by using the optimality conditions for the regions TZl-2, T^ls, • • • , 
respectively, in the reverse direction. These equations can be solved by appropriate nonlinear equation 
solvers and do not require the use of large number of training samples. Thus AQPA is significantly faster 
than GLA and is applicable to the case of large number of feedback bits. Note however, as this is an 
approximate algorithm only, the performance of this algorithm becomes comparable to modified GLA 
only for large values of L. Numerical results presented in the next section illustrate that AQPA performs 
extremely well for L > 16. 

V. Optimum Quantized Power Allocation with Noisy Limited Feedback 

In the previous section, we assumed ideal error-free feedback in the limited feedback model. How- 
ever, feedback channel noise can result in unavoidable erroneous feedback, which can cause the SU- 
TX incorrectly selecting an incorrect transmission strategy and thus dramatically degrade the capacity 
performance. In this section, we allow noise in the limit feedback channel model and study the ergodic 
capacity maximization problem Q with noisy limited feedback. The noisy feedback link, assumed to 
be memoryless, is characterized by the index transition probabilities pkj,(k,j = 1,...,L), which is 
defined as the probability of receiving index k at the SU-TX, given index j was sent from the CR 
service provider/band manager. Given B = log 2 L bits feedback, denote binary representation of index 
k and j as kik-z . . . ks and jxj^ ■ ■ -3b respectively, where k n , j n 6 {0, 1} for n = 1, ... B, and fci, j\ 
represent the most significant bit. We model the noisy feedback channel as B independent uses of a 
binary symmetric channel with crossover probability qf for every feedback bit. Since bit errors are used 
to be independent, p kj = Yln=i Pk n j n — <7/ fe ' J (l — qf) B ~ dkJ , where dkj is the Hamming distance between 
the binary representations of k and j lfT8l|[T9l . 
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Thus problem Q with noisy limited feedback can be reformulated as 

L L AI 

3=1 fc=l i=l 
L I L I ^ Af 

J2^ E ^PMPkjM^j) <QL 9 , Vz, EE^m E E ^l^)^ Pr ^^ F ^ (13) 
j=l fc=l j=l fc=l 1=1 

Note that the binary codewords representing the feedback indices for a power codebook of size L can be 
designed in L\ different ways. Thus finding the optimal index assignment can be done by an exhaustive 
search for small B. For large B, one could resort to some low-complexity suboptimal index assignment 
schemes like lfT9l . Note that such index reassignment schemes will yield the same codebook but with 
its power vectors in different location [19]. Here, given a fixed index assignment scheme, we simply 
concentrate on finding the optimum CSI partitions and power codebook V that jointly optimizes 

the ergodic capacity of SU under the long term average transmit power constraint and average interference 
constraint given by ( [T3~j ). 

Again, to keep things simple, we look at narrowband spectrum- sharing case (M=l). Using the simplified 



notations pj, j = 1, 2, . . . , L, and gi, g , Q avg , we write the Lagrangian for the problem ( 13 ) with M = 1 
as 

L L 

- H90Pk\'R-j]pkjPr(Tl j ) + XP avg + flQavg (14) 

3=1 k=l 

where A and // are the nonnegative Lagrange multipliers associated with the ATP constraint and AIP 
constraint respectively. Thus the Lagrange dual function g(\, fJ>) is defined as 

L L 

^"5^ U .EE E t lo g( 1 + 9iPk) - Xp k - ng p k \'R'j}pkjPr(n j ) (15) 

3=1 k=l 



and the corresponding dual problem is min A >o, At >o5'(A,/i) + XP a v g + pQ 



avg- 



We can solve the optimization problem (15) with fixed A and fi using another modified GLA, (termed 
as modified GLA-2 to distinguish it from the noise free case) by repeating the following two steps 
until convergence: 1) Using large number of training samples for (g ,gi), assign individual (go,gi) 

samples to Kj if 52k=iQ°&0- + 9iPk) ~ ^Pk - P9oPk)pkj > Efc=i( 1 °g( 1 + 9iPk) ~ >>Pk ~ P9oPk)pkn, n = 
1, . . . L, n 7^ j, Vj = 1, . . . , L. 2) Given a partition, the optimal power codebook is given by solving 
the convex optimization problem argmax pfc>0 J2j =1 -E[k)g(l + g±Pk) — Xpk — P9oPk\T^j]PkjP r (J^j), Wk = 
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1,2, ... ,L. One can then obtain the optimal power as max(p£, 0), where p* k is the solution to the equation 

ZU E[u^r k - (A + m)mp* = o- 

For this power codebook, the optimal values A and fx can then be obtained numerically by solving 
the associated KKT conditions. One can repeat the modified GLA-2 and the algorithm for finding A, ii 
iteratively until a satisfactory convergence criterion is met. The extension to the wideband case is obvious 
and is thus omitted. 

VI. Numerical Results 

In this section, we will evaluate the performance of the designed power allocation strategies via numer- 
ical simulations. We implement a wideband spectrum sharing system with one SU and M independent 
frequency bands (each band is originally licensed to a PU), where all the channels involved are assumed 
to undergo Rayleigh fading, namely all g and g 1 are exponentially distributed with unit mean. For each 
simulation, 100,000 randomly generated channel realizations for each g or g 1 are used. 

Fig. [4] shows with prefect CSI, the capacity performance of SU-TX, which shares spectrum with four 
PUs (M=4), with four different AIP constraints thresholds, i.e, (Q av ±, Q a v2, Qav3, Qavi)=(— 5 dB, —5 dB, 
dB, dB), (Qavi j Qav2, Qav3, Qavd<0 dB, dB, dB, dB), (Q avl , Q av2 , Q av3 , Q „ 4 )=(-5 dB, dB, 
dB, 5 dB) and (Q a vi,Qav2, Qav3, QavA)={—^ dB, dB, 5 dB, 5 dB). An interesting observation from Fig. 
[4] is that when P av is small (P av < —5 dB), no matter what the value of (Qavi, Qav2, Qav3, Qavi) is, the 
capacity performance of four curves are almost indistinguishable. This is due to the fact that (see Theorem 
[I]), when P av < mm(Q avl ,Q aV 2, Q a v3, Qavi) (since g 1 is i.i.d), all AIP constraints become inactive. As the 
value of P av increases, the capacity performance with different (Q av \, Q aV 2, Qav3, Qavi) gradually becomes 
distinguishable, since in this case, the ATP and at least one AIP constraint are effective. However, as 
P av increases beyond a certain threshold, the capacity curves start to saturate, due to the fact that when 
Pav > lEtiE[(^-r - where ^ is given by solving i<M^ - ^) + ] = Q av , only the AIP 

constraints are active. Thus no matter how P av changes, if (Q av i, Qav2, Qav3, Qavi) are fixed, the capacity 
will be unchanged. A similar observation for a narrowband spectrum sharing model with full CSI was 
made in [6J. One should note that theoretically, the ATP corresponding to the optimal power allocation 
law maximizing the SU ergodic capacity over a Rayleigh fading channel under an AIP constraint with 
perfect CSI is infinity IfTOll . Since here we use large numbers of randomly generated channel realizations 
samples in the simulation studies, the ATP for maximizing SU ergodic capacity under an AIP constraint 
is large but not infinite. 
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Fig. [5] shows the capacity performance of SU sharing a narrowband spectrum with one PU with limited 
feedback for Q av = —5 dB and Q av = dB respectively, and illustrates the effect of increasing the number 
of feedback bits on the capacity performance. For comparison, we also plot the corresponding capacity 
performance with full CSI. The striking observation from Fig. [5] is that introducing one extra bit of 
feedback substantially reduces the gap with capacity based on perfect CSI. This property is not very 
obvious when P av is small, for example when P av < —5 dB (P av < dB) for Q av = — 5 dB (Q av = 
dB). But with increasing P av , it becomes more pronounced. To be specific, for Q av = —5 dB case, at 
P av = 10 dB, with 1 bit, 2 bits and 3 bits of feedback, the percentage capacity loss is approximately 
21.23%, 6.21% and 1.62% respectively, and for both Q av = —5 dB and Q av = dB cases, only 3 
bits feedback can result in secondary ergodic capacity very close to that with full CSI. This is very 
encouraging since only a small number of bits of feedback are required to achieve close performance to 
the full CSI case. It can be also seen that the capacity performance with large AIP threshold (Q av = 
dB) outperform the ones with low AIP threshold (Q av = —5 dB), as expected. A similar behaviour can 
be also observed in Fig. [6] for a wideband spectrum sharing case (M = 4)((Q avl , Q aV 2, Qav3, Qav&) = (— 10 
dB, -5 dB, dB, 5 dB)). 

In Fig. [7] we compare the performance of AQPA with modified GLA, where SU shares the spectrum 
with four PUs (M = 4) and AIP constraint thresholds (Q av i, Q aV 2, Q a v3, Qavd = (— 10 dB, —5 dB, dB, 
5 dB). It illustrated that with the same number of bits of feedback, the gap between AQPA and modified 
GLA becomes smaller as L increases. For example, when P av = 15 dB, the capacity loss by using AQPA 
instead of GLA is about 8.38%, 3.12% and 1.42% for 2 bits, 3 bits and 4 bits feedback respectively. It 
is clearly seen that AQPA with 4 bits feedback can almost approach the full CSI performance. It is also 
noticed that for a fixed A and /i with M = 4 and 4 bits of feedback, AQPA is approximately 10 times 
faster than GLA operating with 100,000 training samples on a Pentium 3 processor. 

Finally, we investigate SU ergodic capacity performance with noisy limited feedback in Fig. [8j for a 
wideband spectrum sharing case (M = 4 and (Q a vi, Qav2, Qav3, Qavd={— 10 dB, —5 dB, dB, 5 dB)). 
It can be observed that as the feedback becomes less reliable (the crossover probability q$ increases), 
significant capacity performance degradation occurs, especially in high P avg - For example, when P avg = 10 
dB, for 3 (2) bits feedback, a noisy feedback channel with = 0.01 and qj = 0.1 can result in 
approximately 3.843% (4.769%) and 17.394% (18.783%) capacity loss respectively, compared to the noise- 
free case. This clearly illustrates that as the quality of feedback link degrades, the benefit of designing 
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an optimal power codebook diminishes rapidly. 



VII. Conclusions and extensions 

We have derived quantized power allocation algorithms for a wideband spectrum sharing system with 
one secondary user and multiple primary users, each licensed to use a separate frequency band, each 
band modelled as independent block fading channels. The objective has been to maximize the SU ergodic 
capacity under an average sum transmit power constraint and individual average interference constraints 
at the PU receivers. Modified Generalized Lloyd-type algorithms (GLA) have been derived and various 
properties of the quantized power allocation laws have been presented, along with a rigorous convergence 
and consistency proof of the modified GLA based algorithm. By appropriately exploiting the properties 
of the quantized power values for large number of bits of feedback, we have also derived approximate 
quantized power allocation algorithms that perform very close to the modified GLA based algorithms but 
are significantly faster. Finally, we have presented an extension of the modified GLA based quantized 
power allocation algorithm to the case of noisy feedback channels. Future work will include deriving 
expressions for asymptotic (as the number of feedback bits goes to infinity) capacity loss with quantized 
power allocation, consideration of primary interference at the secondary receiver and designing of optimal 
index assignment schemes for quantized power allocation with noisy limited feedback. 

Appendix 

A. Proof of Theorem [7] 

1) Note that the Karush-Kuhn-Tucker (KKT) conditions are necessary and sufficient for a convex 
optimization problem. This implies that all the conditions stated in Theorem [TJ are necessary and sufficient. 
When = 0, from the complementary slackness condition, the constraint i Y^fLi -^[Pi(So> 8i)] — Pavg 
does not come into play. In this case, the optimization problem d2b becomes M completely independent 
parallel subproblems all having the same structure: 

max E[\og(l + gipiig^gj)} 

Pi(g :gl)>° 

s.t. E[gl Pi (g , gl )] < Ql g , Vz = 1, 2, . . . , M (16) 

and it is easy to verify that in the above optimization problem, each constraint holds with equality, namely 
-^[fl , oP*(§0' Si)] = Qavg ^z. Thus for each i, from the complementary slackness condition, one can easily 
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show that ji{ > 0. Hence, in this case, we have the optimal solution 



P*(go,gi) = (— -~) + V* (17) 
K9l 9i 



where /if is determined such that E[g l (—^-— 4) + ] = Q l aV q From feasibility, we also have jr Y^iLi -^[("TT - 



\9l 

,)+]< P-..„. 
Si 



2) When > 0, again from the complementary slackness condition, i YltLi -^[Pi(go> 8i)] = 

• If /if > 0, then corresponding AIP constraint must satisfy with equality (i?[(?Qpj(g , g 1 )] = Q l avg ) 
and hence the optimal solution for the z-th channel is 

K(g ,gi) = (— ^r-- i) + (18) 

where /if is determined from #[#0( 3^^7^ - ^) + ] = given A / . 

• If /if = 0, then the corresponding AIP constraint satisfies £?[5&Pi(g , gj)] < Qa^ g > an d m this case 
the optimal solution for the i-th channel is 

P?(go,gi) = (i-V (19) 

In this scenario we also have E[g l Jxj — 4) + ] = -^[(yj — 4) + ] < Qaw since g and g 1 are independent, 
and E[gl\ = l,Vi. 

Thus when A^ > 0, the optimal solution is given by 



, - 7) + if - ^) + ] < 

P<(go,gi) = < 1 1 < 2() > 

otherwise 



where A-^ is determined such that Yl^=x E\Pi] = Pavg- 



B. Proof of Theorem [2] 

Proof: For the modified GLA, one can define a distortion measure d((g , gi),p) = — (log(l+gip) — 
\p — figop). For such non-difference distortion measures, following [8J, one can ensure nonnegativity of 
the distortion measure by introducing a modified distortion measure as d((go, gi),p) = d((go, 9i),p) — 
mm p d((g , gi),p). Since d((g , gi),p) is a convex function of p for fixed (go,g-\), we get the unique 
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minimum p* = (x^-^) + > thus min P d ((9o, 9i),p) = d{{9o, 9i),P*)- Therefore we have d((g , g x ),p) > 
0. Since d((g , gi),p*) is constant for a given (go,gi), thus using distortion measure d((go, gi),p) instead 
of d((go, gi),p) does not affect the results of modified GLA. One can easily show that d satisfies the 
following properties: (1) d is continuous and d E [0, oo), (2)d((g , g\),p) is a convex function of p for 
each fixed (g ,gi), (3) for each (go, gi),d((g , gi),p) oo, as (g ,gi) -> (g ,gi) and || p ||-»- oo, and 
(4) the partition boundaries in the channel space (go,gi) have zero probability. 

Properties 1), 2) and 3) are easy to show and the proofs here are omitted. Property 4) holds due to the 
assumption of continuous fading channels in this work. Note that this is also a necessary condition for 
a codebook to be optimal for a given partition lfT4"l . Note also that the popular fading distributions such 
as Rayleigh, Rician and Nakagami and Log-normal etc. all satisfy the absolutely continuity assumption. 
It is then easy to show that for these types of fading scenarios, the cumulative distribution function (cdf) 
of (go,gi), denoted by F, satisfies the following properties lfT5ll : (5) F contains no singular-continuous 
part and (6) J d((go, gi),p)dF(g , gi) < oo for each p (implying a finite average distortion). Next, let 
g denote (go,gi). Noting that {g(w)} is a stationary ergodic sequence with a cdf F, and letting F nu) be 
the empirical distribution function of the first n members of the sequence lfT5l . one can show that for 
almost every u, and F satisfy (see Lemma 4 of lfT5l ) (7) {F n } converges weakly to the F and 

(8) lim n / d( (g , g 1 ),p)dF n (g ,gi) = J d((g Q , g 1 ),p)dF(g , g t ), for every p. 

Hence, from [fT51 . we can conclude that the modified GLA satisfies properties 1) to 8). Therefore, 
Lemmas 1-3 of [15! are applicable to the modified GLA with probability one and the modified GLA 
satisfies the global convergence and empirical consistency properties as defined in |[T5i 

C. Proof of Theorem [?] i) 

Proof: We need to prove that for any two adjacent regions IZj and lZj + i,j = 1, . . . ,L — 1, 

log(^2_) 

Pj > Pj+i- Given an arbitrary g satisfying < g < ~( p — — ^) (assuming ji > 0), suppose there 
is a point (go,gf) E TZj and a point (<7o,<7i) E TZj+i (neither of these two points is on the boundary), 
and let (go,g\) denote the point on the boundary corresponding to the same go, which from Lemma [T| 
is given by g\ = * e v+ l >g Hp j - P]+1 ) Then ' we have 9i > 9i > 9v Now suppose pj < p j+1 . Since 
(9o, 9i) e TZj, we have log(l + g\pf) - Xpj - fJ,g Pj > log(l + gftpj+i) - Xp j+1 - ngoPj+i As pj < p j+1 , 
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we have 



(A + A*0>)(Pi+i-Pi) > ^ \ + fl +1 ) 

1 "I" yiPj 

e (.*+Wo)(Pi+i-Pi) _ 1 > ^J(p i+1 _^. e ( A +M9o)(Pi+i-Pi)) (21) 

We akn have n b - e (A+ " g o )(p J- p J+i ) -l _ e^+^oX^+i-^)-! , irrmlies -(A+Wo)(w+i-w)_ 

we also nave ^ - p ._ pj+ie (\+ M0 )(p j - Pj+1 ) - p j+1 -p ]e ^+^ )(P J+1 - Pj ) • iNote mat Pj+i > Pj implies 



1 > 0. Since g\ > 0, we have pj + i — p J e < ^ A+M5,0 ^ Pj+1 p ^ > 0. Applying the above result to (21 ), we obtain, 
< — - - — (A+ Mg0 )( Pj - + 7-pj) — <7i which is a contradiction to gf > g°. Similarly, we can also prove that if 
Pj < Pj+i, we have g\ > g\ which is a contradiction to g\ < g\. Thus we must have pj > Pj+i. 

D. Proof for Theorem [3] ii) 

Proof: From Lemma 2, the boundary between any two adjacent regions IZj and TZj+i is given by 

e (X+fj,go)(pj-p J+1 ) _ -y e (\+fj,g )pj _ e (X+figo)p j+ i 

9i 



Pj — pj +ie (\+/J,go)(pj-Pj+i) p^ e (X+ngo)p J+ i — pj +ie (X+fig )pj 
JX+/ig )p € ( _ \ 

= (A + M^o) ,1 \ n! \ >* + m (22) 

where the last equality follows from the mean value theorem for some p € E (pj + i,pj). The last inequality 
holds since we have p je ^ x+ ^ 9a ^ > p^+^Pj+i an d - Pj+ie &+Mo)p, > - Pj+ie {>>+w>)Pi m By rearranging, 

we Pet e (>+» 8 o) Pe(Pj -p 3 . +1 ) . I 

E. Proof of Theorem \3\ Hi) 

Proof: Given a fixed channel partitioning scheme, the optimal quantized power for IZj is obtained 
as pj = max(p*, 0), Vj, where p* is determined by solving the equation E[ j^^-_ — (A + fig )\TZj] = 0. 
We can see that if E[gi\1Zj] < E[(X + fig )\1Zj], then to satisfy the equation, p* < 0, implying pj = 
max(p*,0) = 0. On the other hand, if E[gi\1Zj] > E[(X + fig )\7lj}, p* has to be strictly positive in 
order to satisfy the optimality equation, implying max(p*,0) = p*. We know from Lemma [3] ii) that 
all boundaries between any two adjacent regions have a lower bound given by g 1 > A + fxg , i.e. for 
any given (go,gi) belonging to any of the first L — 1 regions, g\ > A + /xpo- Thus for the first L — 1 
regions, E[gi\7Zj]Pr{7Zj} > E[(X + fig )\1Zj}Pr{7lj} Therefore the optimal quantized power in the first 
L — 1 regions is strictly positive. This cannot be said however for p L as for 1Z L , we cannot guarantee 
gi > A + ng for any given (g , gi) pair in that region. It is thus possible to have p L to be zero. The next 
result shows under what circumstances one can have pl to be exactly 0. 
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F. Proof for Theorem ^iv) 

Proof: 1) We know from Theorem |3j iii) that we always have £'[ 1+ ^ 1 ip — (A + ngo)\lZj] = 0, j = 
1, . . . , L — 1, and for the region TZl, this equation may not satisfied when p L = 0. Let us assume that 
PL > 0. Then we have ££ =1 E l x + ^o|^]Pr{^} = £[i^l^i]M^i}> implying A + pt = 

EU E lTT^-\ n ^ Pr ^} < EU E \9i\^)Pr{^} = 1 since ZU^ilKAPriKj} = E[ 9i ] = 1, 
for i = 0, 1. Hence if A + n > 1, we must have p L = 0. 



From the optimality equation, one can write p,i = x +nE[gl\n] wnen Pi > 0' ^ * s obvious that p, < 



x+nE[go\iz ] ' ^ = 1, 2, . . . , L — 1. Since p L > 0, this is also true for region Therefore when /i / 0, 
A*<U = /iEtiP^hl^]Pr(^) < Ef=i A+SllTA] Pr (^) < ££=i = I- Similarly, if A ^ 0, 

^Pavg < 1- Thus /i > 1 implies Q a „ < 1 and A > 1 implies P av < 1. 

2) Next, we will show that no matter what A, /i is, pi must be zero for a sufficiently large L and 
Hindoo p L _i = . 

(1) First, we will prove that as L — > oo, the boundary between TZl-i and TZl approaches its limiting 
boundary g x = , where 5* = \im L -+ooPL- Given pi > ■ ■ ■ > p L > 0, it is clear that the 
sequence {pi}, i = 1,2, ... ,L is a monotonically decreasing sequence bounded below, therefore 
it must converge to its greatest-lower bound 5* (5* = lirn^oo pl > ) as L — > oo. Therefore, 
it can be easily shown that for an arbitrarily small e > 0, we always can find a sufficiently large 
L such that pi-i — Pl < e- Thus, as L — > oo, (pl-i — Pl) 0. Using this result, we can show 
that the boundary between TZl-i and TZl approaches the limiting boundary gi = 1 _^+^g° )s* ( or 
A + /jgo = x+g lS * ) as L — > oo, (since this boundary can be written as A + jjg = — V ^-^ L — » anc ^ 

log( 1+91?L - 1 ) 

lim i ^^(lim (pi _ 1 _ Pi) ^ ) = lim^ IT g^ = T ^ F ). 

(2) Suppose there exists a (A,//) such that pi > for any arbitrarily large L (implying S* > 0). 
Thus for any L, p L satisfies E[ 1+ 9 ^ — (A + fig )\lZ L ] = 0. From (1), we have as L — > oo, the 
boundary between 1Zl-i and TZl approaches its limit A + p,g = j^^- Note that for a finite 
value of L, the region TZl can be divided into two parts TZli and TZl2 where TZl\ corresponds to 

log( l+SlP£-l ) 

— p^-pl — < A + fig < 1+ 9 g \ 5 , and TZ L 2 corresponds to 1+ 9 g iS » < A + ngo < oo. As L becomes 
arbitrarily large, the region TZli becomes vanishingly small, and one obtains E(\ + ngo\TZL) > 
E ( i+ g \s* > E ( i+g\p L \^ L ) ^ or a sufficiently large L, which is a contradiction to the KKT 
optimality condition for p L > 0. Hence no matter what A, fi is, pl must be zero for a sufficiently 
large L. And 5* = lim^oo = 0. 



23 



(3) 5* = implies the boundary between TZl-i and 1Z L approaches g ± — A + /^go as L — > oo, and 
since as L — > oo, (pl-i — Pl) — >■ and p L = 0, we have Hindoo pl~i = 0. 
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Fig. 1. System model for wideband spectrum-sharing scenario. 
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Fig. 4. SU ergodic capacity for spectrum sharing with 4 PUs and prefect CSI at SU-TX 
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Fig. 5. SU Ergodic capacity with quantized power allocation (using GLA) with one PU for Q a 
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Fig. 6. SU Ergodic capacity with quantized power allocation (GLA) with four PUs (M = 4) 
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Fig. 7. Capacity performance of AQPA with four PUs (M = 4) 
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Fig. 8. Capacity performance of noisy limited feedback with four PUs (M=4) and different BSC crossover probabilities 
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